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ABSTRACT 



Recent studies revealed a strong correlation between the star formation rate (SFR) and stellar mass of the majority of star-forming 
galaxies, the so-called star-forming main sequence. An empirical modeling approach (the 2-SFM framework) which distinguishes 
between the main sequence and rarer starburst galaxies is capable of reproducing most statistical properties of infrared galaxies, such 
^^ . as number counts, luminosity functions, and redshift distributions. In this paper, we extend this approach by establishing a connection 

between stellar mass and halo mass with the technique of abundance matching. Based on a few, simple assumptions and a physically 
motivated formalism, our model successfully predicts the (cross-)power spectra of the cosmic infrared background (CIB), the cross- 
correlation between CIB and cosmic microwave background (CMB) lensing, and the correlation functions of bright, resolved infrared 
O p galaxies measured by Herschel, Planck, ACT and SPT We use this model to infer the redshift distribution of CIB-anisotropies and of 

the CIBxCMB lensing signal, as well as the level of correlation between CIB-anisotropies at different wavelengths. We study the link 
between dark matter halos and star-forming galaxies in the framework of our model. We predict that more than 90% of cosmic star 
formation activity occurs in halos with masses between 10"^ and 10"-^ Mq. Taking into account subsequent mass growth of halos, this 
implies that the majority of stars were initially (at z >3) formed in the progenitors of clusters (M/,(z = 0) > 10'^^^ Mq), then in groups 
(10''-^ < Mhiz = 0) < 10'^^-^ Mo) at 0.5 < z < 3 and finally in Milky-Way-like halos (10"^ < M,,(z = 0) < lO'^-^ Mq) at z < 0.5. At all 
redshifts, the dominant contribution to the star formation rate density stems from halos of mass ~ 10'^ Mq, in which the instantaneous 
star formation efficiency - defined here as the ratio between SFR and baryonic accretion rate - is maximal (~70%). The strong 
redshift-evolution of SFR in the galaxies that dominate the CIB is thus plausibly driven by increased accretion from the cosmic web 
onto halos of this characteristic mass scale. Material (effective spectral energy distributions, differential emissivities of halos, relations 
between Mi, and SFR) associated to this model is available at http://irfu.cea. fr/Sap/Phocea/Page/index.php?id=5 37] 

Key words. Galaxies: star formation - Galaxies: halos - Galaxies: statistics - Cosmology: diffuse radiation - Cosmology: dark 
matter -Submillimeter: galaxies 



'1. Introduction l2004l) . only a small fraction of this background can be directly 

resolved into individual sources at w avelengths longer than 

^ IA detailed understanding of galaxy formation in the cosmologi- 250/zm dOliveretal .1120101: iBethermin'e t al. 201 23). where the 

• fH cal context is one of the main problems of modern astrophysics. CIB becomes dominated by z>l so urces jLagache et al.ll2003l 

/^ The star formation rate (SFR) of galaxies across cosmic time is |2005l; iBethermin et al.ll201 U l2012cl) . We thus have to study the 

H .one of the key observables to understand their evolution. How- statistical properties of the unresolved background, if we want 

. 5t -ever, measurements of SFR are difficult, because most of the UV to unveil the infrared properties of galaxies which host the bulk 

light emitted by young massive stars is absorbed by interstellar of the obscured star formation at high redshift. 
dust. This light is re-radiated in the infrared between 6 /im and 
1 mm. The cosmic infrar ed background (CIB), detected for the 

first time in FIRA S data (iPuget et al.lll996l : iFixsen et al.lll998l: We can use statistical tools to measure the pho tometric 

iHauser et al.lll998l) . is the relic of all dust emissions since the properties of galaxies emitting the CIB. P(D) analysis dCondonl 

recombination. It is the strongest background after the cosmic 1974; Patanchonet al. 2009) is a method measuring the flux 

microwave background (CM B) and contains half of the energy distribution of sources below the confusion limit by consid- 

emitte d after recombination dHauser & Dwekir200lt iDole et all ering only the pix el histogram of an infrared/millimeter map. 

|2006|) . Identifying the sources responsible for this background iGlenn et al] ( 1201 Ol) managed to measure the number counts (flux 

and their physical properties is thus crucial to understanding distribution) of SPIRE sources down to one order of magnitude 

the star formation history in the Universe. Unfortunately, below the confusion limit using P(D) and resolve ~2/3 of the 

because of the confusion caused by t he limited resolution of the CIB into individual sources. Strong constraints on co ntributions 

cuiTent infrared/millimeter facilities (ICondonlll974t iDole et all to the CIB were also derived by stacking analyses (iDole et al.l 
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l2006t iMarsden et al.l l2009l) . This method allows to measure 
the mean flux of a population individually detected at a shorter 
wavelengtlubut not in the far-infrared/millimeter, by stacking of 
cutout images centered on short-wavelength detections. Number 
counts below the c onfusion limit were me asured with a method 
based on stacking (iBethermin et al.l2010alfa) . In addition to this, 
iBethermin et alJ (1201 2cl) also measured counts per redshift slice 
using an input catalog containing both 24 /im fluxes and redshifts 
and a complex reconstruction of the counts based on stacking. 
These analyses provided constraints on the CIB redshift dis- 
tribution. Nevertheless, an empirical model is still necessary 
to deduce the obscured star formatio n history from number 
counts at various wavelengths (e.g. JLe Borgne et al. 20091: 



Vahante et al. 20091: Franceschiniet al. KOlOt Bethermin et al.l 



20111: IMarsden et al.1 l20lUTRahmati&vanderWerj 120111: 



Lapi et al.ll201 iblGruppioni et al.ll201 H) 



Large sca l e CIB anisotropics me a sured by Spitzer 
(iLagache et al.l 120071 iGrossan & SmootI 120071: iPenin et al.1 
l2012bt) . BL AST (IViero et al.l l2009l) . "flie South Pole Tele- 
scope (SPT,lHalletal J 120101) Herschel (JAmblardet al. 2011, 
IViero et al.ll2012h and Planck (iPlanck Collaboration et aDl201 if 
also provide degenerate constraints on the evolution of infrared- 
galaxy emissivities and the link between infrared galaxies and 
dark matter halos. This degeneracy can be broken combining 
anisotropy information with infrared number counts (see above). 
The first generation of models used to predict/interpret CIB 
anisotropics was based on a combination of an evolutionary 
model of emissivities of infrared galaxies and a linear bias or 
a halo occupation distributio n (HOD) model describing the 
spatia l distribution of galaxies (Knox et al. 20011: Lagache et al. 



2007|: lAmblard & Cooravl l2007t IViero et al.l l2009t iHal et al 



20 lot IPlanck Collab oration et al.1 120111: lAmblard etal J 120111: 



Penin et aP l2012at ,Xia et ^ l2012h . The emissivities 



deduced from a model of galaxy evolution (e.g. IBethermin et aP 
1201 U iLapi et al] 20111) or r epresented by a simp le parametric 



function (iHall et al.ll2010t lAmblard et all l201ll) . However, 



these models assume that there is no dependency between 
clustering and luminosity and in general a single HOD or 
linear bias for all redshifts. Consequently, these models have 
difficulties to fit all wavelengths simultaneously. IShang et al] 
(1201 2|) proposed a new approach assuming an infrared-light- 
to-mass ratio which varies with halo mass and redshift (see 
also the De Bernardis & Coorav 2012 approach focused only 
on 250jum). This new model is also able to roughly reproduce 
the number counts (LFs respectively), though their descrip- 
tion of infrared galaxies is simplistic (a single SED for all 
galaxies at all redshifts, no scatter on the mass-to -light ratio). 
Another approach was proposed by lAddison et al.l (12012,) . who 
co mbine a backwar d -evolu tion counts model very similar to 
the IBethermin et al.l ( 1201 lb approach and a scale-dependent 
effective bias of infrared galaxies to predict the CIB power 
spectrum. This simplified approach is very efficient in fitting the 
data, but is purely descriptive and provides little information on 
the physical link between galaxies and dark matter halos. 

In this paper, we propose a new approach to model both CIB 
anisotropics and galaxy number counts based on the observed 
relation between physical properties in galaxies and their evolu- 
tion with redshift. We use the stellar mass (M*) as a gateway 



' 24fim is often used because ~ 80% of the background is resolved 
into so urces at this wavelength jPapovich et alj[2004l : IBethermin et al.l 
I2010ah . 



to link the halo mass (M/,) and star formation rate (SFR). The 
stochastic link between SFR and M* is modeled following the 
Bethermin et al. (2012a) model (hereafter B12), which is based 
on the observation of a main-sequence (MS) of star-forming 
galaxies (i.e. a strong correlation between stel lar mass and 



star f o rmation rate SFR evolving with reds h ift. iNoeske et al 



2007 



; lElbaz et al.l l2007t iDaddi et al.l [2007t iRodighiero et~ 



2011|) and provides one of the best fits of mid-infrared-to-radio 
number counts. It contains two types of star-forming galaxies 
with different spectral energy distributions (SEDs): secularly 
star-forming MS galaxies and a population of episodic, probably 
merger-driven starbursts with a strong excess of SFR compared 
to the main sequence following t he 2 star formation modes 
(2SFM) formalism introduced in ISargent et al] (|2012|) . The 
relation bet ween stellar and halo mass is derived by abundance 
matching (IVale & OstrikeJ 12004 IConrov & Wechsled 120091: 
iBehroozi et al.l 1201 Ol) assuming a monotonic relation w ithout 
scatter between these quantities. IBethermin et al.l ( 1201 2bl) used 
this technique to conne ct SFR and M;, for MS galaxies in a 
qualitative way (see alsolWang et al.l2()12b. This pape r improves 
and extends the approach of Bethe rmin et al.l (1201 2bl) in order to 
predict the anisotropy of the CIB and the clustering of infrared 
galaxies. This new formaUsm also permits us to describe in 
a phenomenological way the quenching of star-formation in 
massive galaxies and their satellites and involves a refined 
treatment of sub-halos. 

In Sect.|2] we present the philosophy of our approach. In 
Sect.[3] we detail the ideas on which our model is based, and es- 
pecially how we assign infrared properties to galaxies hosted by 
a dark matter halo. In Sect.|4]and|5] we describe the formalism 
used to compute the power-spectrum of CIB anisotropics and 
angular correlation functions of infrared galaxies, respectively. 
In Sect.|6] we present the results of our modeling and an 
extensive comparison with observations. In Sect|7] we discuss 
the successes, but also the limitation of our model. In Sect.jSJ 
we detail the properties of the CIB predicted a posteriori by 
our model, such as redshift distribution or correlation between 
bands. In Sect.|9] we discuss how the history of star formation 
history proceeds depending on the mass of dark matter halos. 
We conclude in Sect.[TOl In AppendixlAl we provide tables of 
conversion from multipole { to angle and from frequencies to 
wavelengths, as we use both conventions interchangeably in the 
paper 

In this paper, we ass ume a WMAP -7 cosmology 
( iLarson et al.1 1201 Ol) and a IChabried ( l2003h initial mass 
function. 



2. Philosophy of our approach 

The majority of previous CIB models were purely phenomeno- 
logical and describe the emissivity and clustering of infrared 
galaxies using a large set of free parameters. This approach 
is useful to derive quantities such as, e.g. the mass where star 
formation is most efficient. However, it is sometimes hard 
to test the validity of these models as a good fit ca n be easily 
obtained considering their number of free parameters. iKim et aP 
( 1201 2h proposed a physical approach based on a semi-analytical 
galaxy formation model, which unfortunately have substan- 
tial discrepancies with the data. We propose an alternative 
phenomenological approach, which represents an intermediate 
solution between a fully empirical and a physical model. We 
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minimize the number of free parameters and build our analysis 
on the observed relation between physical quantities (e.g. the 
specific star formation rate in main-sequence galaxies calibrated 
from optical, near-infrared, far-infrared and radio data). We thus 
do not aim to fine tune the various parameters of the model, 
but to test if the scaling-laws measured from external datasets 
(measurements on small samples extrapolated to all galaxies, 
optical/near-infrared measurements of stellar mass or SFR) are 
compatible with the data under different scenarios. For this 
reason, we choose an approach based on abundance matching 
with no free parameters to describe the link between stellar 
mass and halo mass and use the B12 model, which follows the 
same philosophy to Unk infrared properties and stellar mass. 



Here Al/, may stand for either the total mass of a halo (Mi,) for 
a main structure or the mass of the sub-halo for a sub-structure 
('Wsub)- This function will be useful for our abundance matching 
procedure presented Sect l3.3l because we will assume that the 
properties of galaxies are linked with Mi,. However, this is not 
exactly a mass function, because sub-halos are counted twice, 
namely both in the total mass function and in the sub-halo mass 
function. Fig.[r| (left-hand column) shows the contribution of 
main and sub-halos at z~0.5. The majority of high-mass halos 
are main-halos, while a large fraction of low-mass halos are 
sub-structures of more massive halos. At A1;,=1O'^M0, ~l/3 of 
halos are sub-halos of more massive halos. 



3. Connecting star formation and halo mass by 
abundance matching 

In this section, we describe how we stochastically assign prop- 
erties of star formation to galaxies as a function of the host halo 
mass, combining prescriptions from several modeling works. In 
Sect. l3.1l and l3.2l we present the halo and stellar mass functions 
used in this paper. We then describe how we connect stellar 
and halo mass by abundance matching in Sect. l3.3l Sect. l3.4l de- 
scribes how we stochastically attribute a star-formation rate from 
the stellar mass using recipes based on the B12 model. Finally, 
we describe how we deduce infrared properties of the galaxies 
from their physical properties in Sect. l3.5l using the B12 model. 

3.1. Halo mass function 

In our analysis, we used the halo mass function (HMF) of 



[Tinker et all ( l2008h (in our notation ^ipf^wy )- We choose halo 
mass (Ml,) to be defined by an overdensity of A = 200 (often 
called M2oo)- This HMF evolves with redshift and was calibrated 
on N-body simulations. We also need the mass function of sub- 
halos, which are supposed to host satellite galaxies. The one we 
adopt here comes from Tinker & Wetzel (2010), and provides 
the mass distribution of sub-halos in a parent halo of total mass 
Mh: 



dN 



0.30 X 



Jlog(msub)xln(10) 

-0.7 

X exp 






(msubIM/,) = 

2.5\ 



-'(^) 



(1) 



where nisuh is the sub-halo mass. In our analysis we neglect sub- 
structures inside sub-halos. The mass function of sub-halos is 
thus 



d^N 



■M,, 



dN 



(ilog(msub) 



(msub = Mh\Mi,) X 



dlog{m^ub)dV 
d^N 



dlog{Mi,)dV 



(Wsub) = 

dlogMh, 



(2) 



We also introduce a pseudo "total" mass function, given by the 
sum of the mass functions of halos and sub-halos; 



d^N 



d\ogiMh)dV 



(Mh) 



d^N 



d\og{M,,)dV 
d^N 



(M,, = Mh) 



dlog(m,ub)dV 



(msub = Mh). 



(3) 



3.2. Stellar mass function 

We used the same stellar mass function (SMF) of star-forming 
galaxies as B12, in order to be consistent with this model which 
is used to link the stellar mass to the infrared properties (see 
Sect. 13.4b . This mass function is parametrized by an evolving 



ISchechteiJ(ll976h function: 



^A^ 



dlogiM^) 



<Pb(z) X 



Ml, 



X exp 



Ml. 

'Ml, 



X ^ln(lO), 
Ml, 

(4) 



with a redshift-invariant characteristic mass Mb and faint-end 
slope umf. <Pb, the characteristic density, is constant between 
z=0 and z=l but decreases at z>l as 



log(<^/,) = \og{(pb){z < 1) + ysFMF(l - z). 



(5) 



The various parameters were chosen to reproduce the observed 
evolution of the mass function of star-forming galaxies. Their 
values are given in B12. 

In order to correctly populate dark matter halos, we also 
need to account for the population of non-star-forming galaxies 
(called quenched galaxies hereafter), which are essentially red, 
passively-evolving, elliptical galaxies below the main sequence. 
Star formation activity in these objects is weak and this pop- 
ulation was thus ignored in the B12 model because of their 
negligible infrared emission. These galaxies howeve r contribute 
signifi cantly to the mass function at high mass (e.g. Illbert et aP 
1201 3l see Fig.[I] upper right panel) and thus are generally the 
kind of galaxy that is encountered in massi ve halos. W e used the 
mass fu nction of qu i escen t galaxies from Illbert et al.l (1201 Oh a t 
z<2 and lllbertet"al] (120131) at z>2, fitted bv a lSchechteii(ll976h 
function. \og{(pi,), log(Mi), and a are interpolated between the 
center of each redshift bins and extrapolated at z > 3. The total 
mass function (showed Fig. [1] upper right panel) is the sum of 
the contribution of both quenched and star-forming galaxies. 
The fraction of quenched galaxies at a given stellar mass and 
redshift is called q'(M^,,^). The upper right panel of Fig.[T]shows 
the mass function and its decomposition into quenched and 
star-forming galaxies. The lower right panel shows the fraction 
of quenched galaxies as a function of stellar mass. At high 
mass (Mi, > 10" Mq) and low redshift (z< 1), the majority of 
galaxies are quenched, when the other regimes are dominated 
by star-forming galaxies. 
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Fig. 1. Illustration of the method used to connect various quantities by abundance matching as described in Sect. l3.3l The dashed line illustrates 
the connection between various quantities. We arbitrarily chose to plot results at z = 0.5. Upper left panel: integral HMF and contribution of main 
halos (dotted orange line) and sub-halos (dot-dash purple line). Upper right panel: integral SMF and contribution of star-forming (long-dashed 
blue line) and quenched (three-dot-dash red line). Lower left panel: variation with halo mass of the fraction of halos that are sub-halos. Lower 
right panel: fraction of quenched galaxies as a function of stellar mass. 



3.3. Connecting stellar mass and halo mass through 
abundance matching 



The abundance matching technique is based on the hypothesis 
of a monotonic link between two quantites. This is a fair as- 
sumption for the link between th e stellar and halo mass of a cen- 
tral galaxy at a given redshift (iMore et al.l l2009t iMoster et al.l 
2010; Behroozi et al. 2010 ). In this work, we will follow 
Behroozi et al.l(l2010h and lWatson & Conrovl(l2013h and assume 
that sub-halos and main-halos follow the same Mi,-Mf, relation. 
Under this assumption, we associate halo mass to a stellar mass 
by putting the n-th most massive galaxy (in term of stellar mass) 
into the n-th most massive halo. In practice, we do not use cat- 
alogues, but analytic mass functions. The non-parametric func- 
tion linking stellar and halo mass (M* = /(At/,)) is thus the 



solution of the implicit equation 

nM,{> f{Mh) ^M^)^viM,„ 



(6) 



where nM,(> M^) is the number density (in comoving units) of 
galaxies more massive than M* (i.e. the integral of the mass 
function), and iimi, the equivalent for halo mass. We can thus 
associate halo mass to stellar mass by taking the halo mass at 
which the number density of galaxies and halos are the same, as 
illustrated by Fig.[T](upper panels). 

In this approach, we neglect the effect of the scatter on 
the stellar-to-halo-mass ratio, which would further complicate 
our analysis of CIB anisotropics. This induces a bias on 
the estimation of /. However, this effect is s maller than the 
statist ical uncertainties for Mi, > 10'"^^ Mq dBehroozi et al.l 
1201 oh . which host mainly passive galaxies (see Fig.[Tl). The 
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scatter around / could also induce a bias on the estimate of the 
observables as the CIB anisotropics. However, the large-scale 
anisotropics are sensitive to the mean emissivity of galaxies and 
are not affected by the scatter. The small-scale Poisson term can 
be computed directly from counts model without assumptions 
on the dark matter (see Sect.H]). The angular correlation of 
bright resolved galaxies can be affected by the scatter on halo 
mass, but the scatter between stellar and halo mass (~ 0.15 dex) 
has the same order of magnitude as the scatter between stellar 
mass and star formation rate, which is taken into account by 
our model (~0. 15-0.2 dex). The impact of the scatter of the 
stellar mass-halo mass relation on the correlation function is 
thus expected to be relatively small. 



3.4. Connecting star formation rate to stellar mass 

In the previous section, we explained how we can assign a stel- 
lar mass to a galaxy knowing its halo mass. Unfortunately, we 
cannot link star formation rate to stellar mass by abundance 
matching assuming a monotonic relation. This hypothesis is 
only valid for main-sequence galaxies, but not for quenched, for 
which sSFR<<sSFRms, where sSFR is the specific star forma- 
tion rate i.e. SFR/M* and sSFRms the typical value of this pa- 
rameter in main-sequence galaxies, and starburst galaxies, for 
which sSFR>>sSFRms- For quenched galaxies, we neglect the 
star-formation and thus take SFR=OM0.yr"' for simplicity. For 
star-forming galaxies (main-sequence and starburst), we assume 
that SFR follows a d ouble log-normal di stribution at fixed red- 
shift and stellar mass dSargent et al.ll2012l B12): 



exp - 



p (log(SFR)|M*,z) oc pms + psB 
(log(SFR)-log(sSFR^jsXM,))2^ 



2^^^ 



+rsB X exp 



(log(SFR) - log(sSFRMs x M*) - Bsb))' ) 



2o-; 



SB 



(7) 



The first term describes the sSFR distribution of main-sequence 
galaxies (pms), and the second one that of starbursts (pss)- Star 
formation in starburst is boosted on average by a factor Bs j3. 
(Tms and ct^b are the dispersion around the central values. The 
evolution of the main-sequence parametrized as in B12: 



sSFRms(z, M^) =sSFRms,o x 



M* 



x(l + min (z, Zevo)) 



7ms 



(8) 



The starburst ratio rse (ie. the relative normalization of the two 
log-normal distributions) is also provided by the B12 model: 



rsBiz) = rsB.o x (1 + min (z, Zs b))^'" , where zs^j = 1. 



(9) 



The values of the parameters in our base model (model A) 
are the ones provided in BI4J. We also used a second version 
of the model (model B) for which the high redshift trend was 



^ A detailed discussion of the SFR-enhancements of starbursts and 
their description in the 2-SFM framework is provided in lSargent et al.l 
(20l3). 

^ Parameters provided in B12 are given assuming a ISalpeten l ll955h 
IMF when this paper assumes a lChabriea l l2003h IMF. A correction of 
0.24 dex have thus to be applied to some of the parameters to take into 
account this difference of IMF. 
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Fig. 2. Upper panel: Evolution of specific star formation rate in main- 
sequence galaxies of 5 x 10' Mq. Lower panel: Evolution of the charac- 
teristic mass (pt, of the mass function with redshift. Both panel: model 
A (B12 model) is represented with a solid line and model B (modi- 
fied version with higher sSFR and lower density at high redshift) with a 
long-dashed line. 



slightly modified, follo wing the findings of a slowly increasing 
sSFR at high redshif t (IStarket al.l 1201 2t Ide Barros et"aDl2012L 
iGonzalez et al.ll2012l) . These higher values are the consequence 
of an improved modeling of the contribution of nebular line 
emissions to flux in near-infrared broad-band filters. For this 
modified version, we assume an evolution of sSFR at z > 2.5 in 
(1 + z). In order to avoid an overprediction of bright millimeter 
counts, we compensate this increase of sSFR by a quicker 
decrease of the characteristic density of the stellar mass function 
in the same redshift range ((1 + z)""*^ instead of (1 H- z)"^''), thus 
keeping the same number of bright objects. As shown in Fig.|2] 
both these scenarios are compatible with the data, because of the 
large scatter on the measurements. In Sect|6] and |7] we discuss 
which scenario is actually favored by infrared observations. 



3.5. Infrared outputs of galaxies 

In massive galaxies, the bulk of the UV light coming from young 
stars is absorbed by dust and re-emitted in the infrared. We can 
thus assume that the bolometric infrared (8 - 1 000 ^m) lumin os- 
ity LiR is proportional to star formation rate dKennicutll 19981 the 
conversi on factor is K - S FR/Lir = 1 x lO^'^Moyr^^L^' if we 
assume a lChabrierl (l2003l) IMF). In low-mass galaxies, a signif- 
icant part of the UV light escapes from the galaxy and infrared 
emission is no longer proportional to SFR. Star formation rates 
can then be estimated from an uncoiTccted UV and an infrared 



Article number, page 5 ofl24l 



component {S FR - S FRuv + S FRir). The infrared luminosity 
is then given by: 



LiR — 



SFR. 



K 



SFR 



-|^Q0.4xruv 



SFR 



xg{M^), (10) 



ruv = 2. Slog 



\SFR 



uv 



■■ 4.07 X log 






39.32. 



(11) 



g(Mi,) tends to at low mass and 1 at high mass. The UV light 
from young stars thus totally escapes the low-mass galaxies, 
but is fully reprocessed and emitted in as infrared radiation 
in massive galaxies. This is due to a larger amount of dust in 
massive galaxies, which causes a larger attenuation. 

We used the same spectral energy distr ibution (SED) tem- 
plates as B12, based on lMagdis et alJd2012l) . There are different 
templates for main-sequence and starburst galaxies, both of 
which are assumed to evolve with redshift. We do not adopt a 
single template for a given type of galaxy at a given redshift, 
but assume a scatter on the mean interstellar radiation field 
(L'^), following B12. These SEDs were calibrated using z < 2 
data, which show a rise of the (rest-frame) dust temperature 
with redshift. At higher redshift, we assume no evolution. 
This assumption will be discussed in Sect.|2] In this paper, we 
neglect the contribution of AGN, which is only significant in the 
mid-infrared (B12). 



4. Computation of CIB power spectrum 

We aim to compare the CIB anisotropics predicted by our 
model to observations in order to test its validity. This section 
presents the formalism used to derive the power-spectrum 
(cross-spectrum) of the CIB at a given waveband (between two 
wavebands, respectively). One of the key benefits of the relation 
we established between SFR and At/, is the fact that we can 
then rely on the well-known clustering properties of dark matter 
halos to predict the clustering of star-forming galax ies, and thus 
of CIB anisotropics. We use a method similar to IShang et alj 
(I2OI2I) . However, we modified their formalism to obtain a more 
natural notation and avoid renormalization of all terms by the 
total emissivity of infrared galaxies at a given redshift. 



4. 1 . Mean infrared emissivities of dark matter halos 

One of the key ingredients to compute is the mean emissivities 
of the halos. Classical CIB models assume that clustering and 
emissivity are independent and compute the total emissivity of 
galaxies at a given redshift. This approximation is not exact 
and both emissivity and clustering vary with halo mass (see e.g. 
iBethermin et al.l2012bl) . We thus introduce the differential emis- 
sivity djv/dlog(Mi,) of dark matter halos as a function of halo 
mass. This differential emissivity is the sum of the contribution 
of central galaxy and satellite galaxies: 



djy 



dlog(M,) 



(Mh,z) 



djv,c 



dlogiMf,) 



(Mh,z) + 



"7v,sub 

dlogiMh) 



(Mi„z). (12) 



1.00 



■o 

0) 



where ruv is the ratio between obscure d and unobscured star ~ 

formation. We use the ruv-M* relation of iPannella et al.l (l2009h E 

to compute §(M*): ° 

I SFR/R 



0.10 r 






0.01 




10 100 

Rest-frame wavelength [wm] 



1000 



Jt 



CO 




10 100 1000 

Observed wavelength [am] 

Fig. 3. Effective SEDs of infrared galaxies used in our model at var- 
ious redshift. The upper panel shows the SEDs in vLy units normalized 
at LiR = 1 Lo as a function of rest-frame wavelength. The lower panel 
shows the ratio between the observed flux density and Lir as a function 
of observed wavelength. 



The contribution of central galaxies to the differential emis- 
sivity is computed from mean infrared flux of galaxies hosted by 
a halo of mass M/, and the HMF: 



djv.c 



dlogiMh) 
SFRms(M 



(M,,z) = 



d^N 



X- 



dlogiMh)dV 
f{Mh),z) 



xD;(l-Hz) 



^ xg{M^^f{Mh),z) 

xsf{z)x(l-q(M^^f(Mh),z)), 



(13) 



where Dc is the comoving distance and sl^(z) is the effective 
SED of infrared galaxies at a given redshift, i.e. the mean flux 
density received from a population of star-forming galaxies with 
a mass corresponding to a mean infrared luminosity of 1 L© 
(see Fig.|3]l, and q (M* - f (Mi,) , z) is the fraction of quenched 
galaxies (i.e. non star-forming galaxies) as defined in Sect. 13. 2! 
The shape of sSFR distribution being invariant with mass in our 
model, the mean flux density coming from a more massive pop- 
ulation can thus be obtained just rescaling this SED taking into 
account the attenuation (^(M*) defined in Eq.fTOt and the SFR- 
M* relation. The effective SED is thus the mean of each type 
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Central + Satellite 
Central only 




halos, the mass function of sub-halos in main halos of this mass, 
and the mean flux density of galaxies hosted by sub-halos: 



Fig. 4. Differential emissivities at 100//m (3000 GHz) and 1.38 mm 
(217 GHz) as a function of halo mass at various redshift predicted by 
model C. The solid lines are the contributions of all the galaxies to the 
infrared flux of a main halo, when the dotted lines indicate only the 
contribution of central galaxies. 



^Jv,sub 



iMh,z)^ 



d^N 



I 



sf(z)- 



d\og(Mi,) 

dN 

-(msubIM/,) X 



X D CI + 7^ 
d\og(Mh)dV '^ ~' 

SFRMs(M„^/(m,ub),z) 

'(ilog(m.sub) """"''"' K 

xg(M^ = /(msub), z) X (1 - ^ (M* = / (m^ub) , z)) dlog(m,^b)- 

(16) 

In this formula, we assume that the quenching of satellite galax- 
ies depends only on their stellar mass (or sub-halo mass). We 
also propose an alternative scenario, where quenching depends 
on the mass of the main halo and not on the mass of the sub-halo. 
In this scenario, satellite galaxies become quenched at the same 
time as the central galaxy in the same parent halo. In practice, 
we replace wzsub by M/, in the last factor of Eq.[T6] which can 
then be moved outside the integral. This scenario is motivated 
by the fact that the fraction of quenched sate lUte gala xies tends 
to be higher in dense environments (e.g. iPark et a l. 2007). 
This phenomenon is often called environmental quenching. 
The modified version of model B where this modification was 
performed is called model C. 

In our computation, the flux densities are not the monochro- 
matic flux densities at the center of the passband filters of each 
instrument, but are computed taking into account the real filter 
profiles. Fig.|4]illustrates the variation of differential emissivities 
with halo mass, redshift and wavelength. Due to the shape of 
the SEDs, long wavelengths have stronger emissivities at high 
redshift. Note that the halo mass dominating the emissivities 
is always ~ 10'^ Mq in agreement with previous works (e.g. 



Conrov & Wechslej 20091: Befliermin et al.ll2012bl: IWang et al.l 
20I2l:lBehroozi et al.ll2012al[ 



of SED weighted by their contributions to the background (pro- 
vided by Eq.|7]l: 



sf( 



z) = r [pus X (if ^) + psB X (sT)] «?sSFR, 



(14) 



where 



<^i 



MS or SB 



/ 



PMSoi SB ({U}\z)xs 



MS or SB 



{{U),z)xd{U) (15) 



with (U) is the intensity of the radiation field (strongly linked 
to dust temperature), p{{U}\z) its probability distribution (this 
encodes the information on the scatter of dust temperatures), 
and iJ^SorSB jj^g ^^^ density received from a Lir = I Lq 
main-sequence or starburst source with a radiation field ([/) 
at redshift z. (iJ:is°'SB^ j^ jjj^^ jjjg average SED of a MS 
or SB galaxies at a given redshift taking into account the 
scatter on (U). Note that the infrared luminosity of these 
effective templates are slightly larger than 1 Lq because of the 
asymmetry of sSFR distribution caused by starburst. The ef- 
fective SEDs from Eq.[T4]used in this work are available onlinqj- 

The contribution of satellite galaxies to differential emissiv- 
ity is the sum of the contribution of all galaxies in sub-halos of a 
main halo of mass M^. It depends on the mass function of main 



|http : //irfu . cea . fr/Sap/Phocea/Page/index . php?id=537 , 



4.2. Power spectrum 

The CIB power spectrum is the 
butions ( Amblard & Coorav 



■ 1200 

Planck Collaboration et all l20 1 iT " 
Penin et alJl2012aHShang et al.ll20I2l) : 



20071: 



su m of three 
Viero et al.l 



contri- 



lAmblard et al. 



2009; 



2011 



- 2-halo term: correlated anisotropics between galaxies in dif- 
ferent halos, which dominates at scales larger than a few ar- 
cminutes. 

- 1-halo term: correlated anisotropics of galaxies inside the 
same main halo, which have a significant impact at scales of 
a few arcminutes. 

- Poisson term: non coiTclated Poisson anisotropics, also 
called shot noise, which dominate at small scales. 



The cross power spectrum of CIB Cfv 
bands (v and v') is thus: 



Cr 



C 



2/1 






between two frequency 
(17) 



The 2-halo and 1-halo terms, which correspond to large 
and intermediate scale anisotropics respectively, are computed 
from the mean emissivities of galaxies and are not affected by 
the stochasticity of the connexions between galaxies and halos. 
The anisotropics caused by the scatter of galaxy properties are 
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described by the Poisson term, which requires no assumption 
regarding the link between dark matter halos and galaxies and 
is deduced only from the B12 model. Each term is calculated 
independently, but based on the same consistent model. 



4.2.1. 2-haloterm 

The 2-halo term is computed from the following formula, sum- 
ming on redshift, but also over all cross-correlation between ha- 
los of various mass: 



The first factor describes geometry. The second one repre- 
sents the various (cross-)correlations between satellite and cen- 
tral galaxies. There is only a factor u for (cross-)correlation be- 
tween satellites and central, because the central is assumed to be 
at the center of the halo and the satellites follow the NFW profile, 
and u^ for satellite-satellite combinations. Finally, we have to 
renormalize by the inverse of the mass function, because the two 
dj/dlog{Mh) factors contain implicitly two ti mes the number 
of ha los, when this should appear only once dCoorav & ShethI 
l2002h. This notation avoid to have to renormalize by jy as in 
i Shangetal.l(l20Tl . 



C 



2/1 



'-ami 

djv',c 



djv.c 



dlog(M'.) 



iz) + 



d\og(Mh) 

^Jv\suh 



(z) + 



"7v.si 



ub 



iilog(M;,) 



(z)m(A:,M,„z) 4.2.3. Poisson term 



dlogiM;^) 
xPUk 



The Poisson term is independent of large-scale halo and only de- 
(z)u(k M' z)\b(Mi 2)fe(M',z)PSiidsonthefluxdistributionof galaxies (number counts). These 



I 



z) d\ogMhd\ogM'i^dz 
(18) 



Poisson ani sotropics for the aut o power spectrum can be com- 
puted from dLagache et al.ll2003l) : 






Jo dSy 



(22) 



The first factor is geometrical. The next two factors contain 
differential emissivities. u{k,Mi,,z) is the Fouri er transform of 
the halo profile assumed her e to be NFW ( iNavar ro et al. 1997). 
In contrast with IShang et aP (1201 2h . this term is placed only in 
front the sub-halo-emissivity term, assuming that sub-halos are 
distributed following the NFW profile. The central galaxy is 
assumed to be at the center of the halo. The clustering term 
(b(Mh)b(M'f^)P\ia] is the cross power-spectrum between halo of 
mass M;, and M'., under the assumpti on that P(k,Mi„M'i^) = 
b{Mh)b{M])P\in dboorav & Shethll2002h . Eq.[T8]can be signifi- 
cantly simplified introducing: 



Here dNjdSy are the differential number counts (see e.g. B12 
for the computational details) and S'v.cut is the flux cut at which 
sources are removed from the maps. The Poisson term of the 
cross-spectrum is slightly more complex to compute: 



^" ^ Jo Jo '^''^^^;^''^'''- 



(23) 



Jv{z,k) 



-^ 



b{M,z) 



djy 



dj. 



,sub 



d\og(Mh) dlogiMh) 



u(k,Mh,z)\dMh, 



(19) 



where d^N/dSydSy' are the multivariate counts (i.e. number of 
sources with a flux between S ,, and Sy + dSyin one band and S J, 
and S'y + dS'y in the other). 

In practice, multivariate counts are hard to compute, and 
summing the contribution of various redshift, type of galaxies, 
infrared luminosity, and radiation field is easier (the derivation 
of this formula is presented in AppendixlBli: 



which is an emissivity weighted by the bias corresponding to 
each halo mass. We can then simplify Eq.fTSl 



C 



2/1 
e.vY' 



r dP, I a 
J dz \D, 



Jyiz, k)J'^iz, k)PUk ^^^z)dz (20) 






"(WU) 



(MS.SB) 



d^N, 



MSo>-SB^2^^^MSorSB(^f^^^^) X s'^^^^'^HU}, z) dL,^d{U} dz 



This way of computing Cf'^^ reduces the number of integrals, 
because Jy can be calculated only once per frequency channel 
and then can be used to derive all the cross-spectra. In addition, 
Jy can also be used to compute the cross-correlation between 
cm and CMB lensing (see Sect.l431l. 



dLiRdV 



-^IR-'v 



(24) 



where L^™^^^({U},z) is the infrared luminosity where the 
source is detected in at least one of the bands (depends on 
redshift, type of galaxy and radiation field), d^Nus or sb /dLlJ^ dV 
is the infrared luminosity function (main-sequence or starburst 
contribution). 



4.2.2. 1 -halo term 

The 1-halo term is computed with: 



"-'■"' ~ J J dz U. 



^Jv,c "7v 



,sub 



dj. 



,sub 



djy 



dlogiMh) d\og(Mh) 



u{k,Mi„z)- 



' dlogiM,) dlogiMh) 

^Jv,suh "Vv'.sub 

d\og{Mh) d\og{Mh) 

X 



d^N 



dlog{Mi,)dV 



4.3. Cross-correlation between CIB and CMB lensing 

In addition to cross-correlations between the CIB in various 
u(k,Mh,z)+ bands, we can also test the predictive power of our model for 
correlation between the CIB and the reconstructed gravitational 
potential derived from distortions of the CMB d ue to gravita- 
tional lensing by large scale halos at z~l-3 (e.g. iHanson et al.l 
|2009[) . This correlation is a direct probe of the link between 
dzdlogMi, the gravitational potential of dark matter halos and infrared 
emission from star-forming galaxies. 
(21) 



u\k,Mh,z)] 
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Because the CMB lensing signal is due to dark matter halo, 
this signal can be modeled with a 2-halo term replacing the 
term of emissivit y by a term linked to the gravita tional poten- 
tial (adapted from lPlanck Collaboration et al.ll2013h : 



C 



t,<!>v 



r dDc I a 
<D(z, () is given by dChallinor & Lewij2005h . 



Uz, miz, k)PUk = — , z)dz, (25) 

Ur 



(26) 



where Hq is the Hubble constant, Om the matter density in units 
of the critical density, and D™^ the comoving distance between 
CMB and us. The impact of sub-halos on CMB lensing (1-halo) 
is quite small and can be neglected for this work. 



5. Computation of angular correlation function of 
resolved infrared galaxies 

In addition to the anisotropics of the faint infrared sources re- 
sponsible for the unresolved background, our model also pro- 
vides predictions of the angular correlations of the individually 
detected bright sources. A formalism taking into account the 
selection function of the resolved sources has to be employed 
and we thus cannot use the same formalism as for the power- 
spectrum. In this paper, we only consider samples selected us- 
ing a flux threshold {S v.cut) at a given wavelength for simplicity. 
In this section, we first explain how we compute the halo occu- 
pation distribution (HOD, i.e. the mean number of central and 
satellite galaxies as a function of halo mass and redshift) for a 
given selection of resolved sources and then how we derive the 
correlation function from the HODs. 

5.1. Halo occupation distribution 

The mean number of central galaxies in a given halo of total 

mass Ml, is: 



{N,)(Mi„ z)^{i-q(M^^f (Mh) , z)) X y r r 

. JsSFR Ja 



type 



<U> 



//(sf s°''S^<t/),z) X LiR (M* = f(M,,) , sSFR) > 5,,eut) 

d{U) dsSFR, (27) 

where q(Mi, - f(Mh),z) is the fraction of quenched objects 
and // is a function having a value of 1 if the condition is true 
and otherwise. This condition is true only if the source is 
sufficiently bright to be detected {Sy > 5 v.cut)- This depends on 
both halo mass, type of galaxy (MS or SB), sSFR, and ([/). 

The number of satellites depends on the total mass of the 
parent halo M;, and is connected to sub-halo via: 



(A'subXM,, 



•Jni 



dN 



,„b '^log('«sub) 



(msublM/,)(A^c>('«sub, z)dm^^b. 



This formula is true only when we assume the same in- 
frared luminosity versus halo mass relation for main and 
sub-halos (model A and B) and consequently the HOD of 
the satellite in a given sub-halo of mass nisub is the same as 



for the central in a main halo of identical mass M/, - m^uh- 
For model C (satellites quenched at the same time as the 
centrals), the factor (A^c)('«sub,z) has to be multiplied by 
(1 - q(M^ = f (Mh) ,z)) / (I - q{M^ ^ f (m^ub) , z)) to take into 
account the fact that in this version of the model the quenching 
is linked to the mass of the central and not the satellite. 

FigH shows the HOD of S 160 > 5 mJy (sources detected 
by the deepest PACS surveys) and ^ggo > 3mJy (typical 
sub-millimeter galaxies) sources at low (z=0.1) and high (z=2) 
redshift. The number of central galaxies decreases very quickly 
at low mass. This sharp cut is due to the selection in infrared flux 
and the rather tight correlation between infrared luminosity and 
halo mass for central galaxies. At low redshift, the probability 
to detect an infrared galaxy in a massive halo is pretty small 
because the vast majority of central galaxies are quenched 
and thus have little infrared emission. The HOD of satellite 
strongly depends on the version of the model. If the quenching 
of satellites is decorrelated from the quenching of the centrals 
(model A and B), the number of detected satellite increases 
quickly with the halo mass. This is not the case for model C, for 
which the satellites are quenched at the same time as the centrals. 



5.2. Angular auto-correiation function 

The angular correlation function (ACF) can be compute from 
the HOD of galaxie s and their redshift distribution (see 
ICoorav & Shethll2002l for a review): 



w(0) 



iXff khP,,^k,z)MkD,0)dz 
[IMdkdzf 



(29) 



where Jq is the zero-th order Bessel function and dN/dz the 
redshift distribution of galaxies, which can be computed from 
counts per redshift slice following iBethermin et al.l (1201 lb . Pgg 
is the sum of two terms corresponding to the clustering of galax- 
ies in two different halos and inside the same halo: 



P,,(k,z)^P^^{k,z) + P'^g(k,z). 



(30) 



Following the standard conventions, we write the 2-halo term as: 



/■2(*.z) 



/ 

J Ml, 



-b(Mh)^ — dMh 



Ml, d\og(Mh)dV «gal 

where (A^gai) = (A^c) + {Nsnb)- The 1-halo term is: 



Pii„(fc,z), (31) 



P'^AKz) 



/ 

Jm 



d^N 2(JVe)(JVsub) + {N,uhfu\k, Mi„z) 

M„ d\og{Mi,)dV «gal 

dMh, 
(32) 



(28) where 



"gal = J 



d^N 



dlog{Mh)dV 



{Ng,i}dMh. 



(33) 
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Fig. 5. Halo occupation distribution, i.e. mean number of detected galaxies in a halo, including its sub-structures, as a function of its total mass 
M;,, of 5 160 > 5mJy (left) and Sgso > 3mJy (right) sources at z = 0. 1 (top) and z=2 (bottom) for central (solid line) and satellite (dashed line) 
galaxies predicted by the model. The HOD of central galaxies is the same for all the version of the models and represented in black. Satellite 
HODs are plotted in blue (red, respectively) for model A/B (C, respectively). 



6. Comparison with the observations 

6.1. CIB power spectrum and number counts 

We compared the predictions of our model with the CIB 
power- spectrum measure d with IRAS (100 u m. 
2012bt) . Spitzer (160/im, iLagache et aP 120071 and iPenin et al. 
2mm . Herschel (250, 350, and 500 um. IVier o et al."2012'). 



Planck (350, 550, 850, and 1 380 um.lPlanck Collaboration et al.i 
120111)). and SPT (1380^m. iHall et al.ll2010h measurements in 
Fig.|6] (left and center). The level of the Poisson anisotropics 
is significantly aff'ected by the flux cut used to mask bright 
sources, and their level increases with the chosen flux cut. The 
correlated anisotropics (1-halo and 2-halo) are less affected by 
the flux cut. When we compare several datasets at the same 
wavelength, we thus compute the model predictions for the 
flux cut of the experiment which is the most sensitive at small 
scales. The flux cuts used to compute the model predictions 
are thus IJy at lOOjum, lOOmJy at 160yum, 0.3 mJy at 250 //m, 
350;um, and 500 //m, 0.54 Jy at 550pm (545 GHz), 0.325 Jy at 
850;um (353 GHz), 20mJy at 1.38 mm (217 GHz). The overall 
agreement with the data is very good for a model not fitted 
on the data. However, some tensions between model A and 
B and data at short and large wavelengths provide interesting 



information. 

The anisotropics at 100/im are dominated by galaxies 
at low redshift (see Sect.[8]l and are thus not affected by the 
evolution of galaxies at z>2.5. Consequently, there is no 
difference between model A and B. On the contrary, model C 
(with quenching of satellite around massive quenched central) 
predicts a lower level of anisotropics at ^ < 5000, which is in 
better agreement with the data, especially around i - 2000. The 
"environmental quenching" in model C reduces strongly the 
1-halo term compared to model B (by a factor of 4 at ^=2000), 
but does not significantly affect the 2-halo term (~15%) and 
does not change the Poisson level at all. However, the former 
difference is not very significant (~ 2cr) and CIB anisotropics 
are obviously not the best probe of such environmental effects. 

The interpretation of Spitze r/MlPS 160jum data is more 
tricky. While the measurements of iLagache et al. (l2007h 
crosses) agree with our data, measurements of iPenin et aP 
2012bl) (asterisks) are a factor 3-10 higher than our model at 
large scales {{ < 2000). Such a large difference is unexpected, 
especially when we consider that the model agrees well with 
the data at 100 and 250 /im. However, MIPS suffers from strong 
1/f noise and is not the best instrument to measure large scale 
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big. 6. Left an d center panels: Clri p ower sp ectrum predicted by our model and comparison with IKAa ( luu um.JHenm et al.lzUl/a . open circles), 
Spitzer (lOOom. lLagache et al.120071 crosses.'Psnin et al.'2012bl. asterisks), Herschel (250, 350, and 500om.'Vi ero et ani2012l . triangles), Planck 



Fig. 6. Left an d center panels: CIB p ower sp ectrum predicted by our model and comparison with IRAS (100 um.|Penin et al.l2012bl . open circles), 

w 

(350, 550, 850, and 1380 um. IPlanck Collaboration et al. (2011), squares), and SPT (1380/jm, Hall et al. 2010, diamonds) measurements. The 
dotted, dot-dashed, and dotted lines represent the Poisson, 1-halo, and 2-halo terms. Right panel: Number counts of infrared galaxies. Data are 
taken from the compilation of measurements in iBethermin et al.l i201lh and B12. Model A, B, and C are represented with a green long-dash, a 
blue short-dash, a red solid line, respectively. 



anisotropies. large scale by a factor 2 at 1.4 mm, while models B and C 

agree with the data, except for a 2cr tension at 1.4 mm in the 

two lowest-multipole bins. However, measuring the CIB at 

Between 250 and 550yum, our model nicely agrees with 217 GHz is difficult because it strongly relies on the correct 

Planck and Herschel measurements. At longer wavelengths subtraction of the CMB. Future analyses will either confirm 

(850yum and 1.4mm), model A over-predicts anisotropies at 
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Fig. 7. Cross-power spectrum between SPIRE bands measured bv lViero et alj ( 12013) and comparison with our model. Model A, B, and C are 
represented with a green long-dash, a blue short-dash, a red solid line, respectively. 
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Fig. 8. Cross-power spectrum between BLAST and ACT measured bv lHaiian et al.l i2011h and comparison with our model. Model A, B, and C 
are represented with a green long-dash, a blue short-dash, a red solid line, respectively. 
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Fig. 9. Cross-power-spectnim between C IB and CMB lensing mea- 
sured by iPlanck Collaboration et alj d2013f) and comparison with our 
model. Model A, B, and C are represented with a green long-dash, a 
blue short-dash, a red solid line, respectively. 
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Fig. 10. Auto-correlation function of various flux-selected sample 
of infrared galaxies and comparison with our model (same color coding 
as in Fig.|6ll. Th e data come from M aglio cchetti et al. (2011) at lOO^m 
and 160 //m and lCoorav eTal] ( 120101) and lMaddo x et al. (2010) at 250, 
350, a nd 5 00 /^m. We used the same flux cuts as Magliocchetti et al.l 
j2011[) and ICoorav et all ( 120101) . Model A, B, and C are represented 
with a green long-dash, a blue short-dash, a red solid line, respectively. 



or refute the presence of this discrepancy. The three models 
agree with the number counts at 850 yum and 1.1mm. The 
difference between model A and B/C is the sSFR at z>2.5 and 
the characteristic density of the mass function of star-forming 
galaxies. Model A assumes a flat sSFR at z>2.5, while models B 
and C assume rising sSFR at z>2.5 compensated by a decrease 
of the characteristic density of the mass function to preserve the 



agreement with the number counts. In these two versions, the 
number of bright objects is thus the same, but the number of 
faint objects is different (see e.g. number counts around 1 mJy at 
1.1 mm). Model A over-predicts the CIB anisotropies because 
of a too high emissivity of high-redshift faint galaxies. 

In the case of the model C (which includes satellite- 
quenching), the 1-halo term never dominates the anisotropies 
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regardless of wavelength or scale. Concerning models A and 
B (without satellite-quenching), the 1-halo term contributes to 
the CIB anisotropics by a similar amount as do the 2-halo and 
Poisson terms at f ~ 2000 at 100 yum, 160 yum, and 250yum. At 
other scales and wavelengths, it represents a minor contribu- 
tion. This result is in contradiction with models assuming no 
mass-luminosity dependence (e.g. iPlanck Collaboration et alJ 
|2011; Penin et al. 2012ai), b ut consistent with the approaches 
assuming a dependence (e.g. lShang et alj|2012h . This reduction 
of the 1-halo term is caused by the fact that satellites are hosted 
by low mass sub-halos and thus have low infrared luminosities. 



6.2. CIB cross-power-spectrum 

In addition to auto power spectra, we can compare the predic- 
tions of our model with cross power spectra. These cross-power 
specta are useful to verify if the level of correlation between 
bands is correct. This also is an indirect test of whether 
SED libraries and redshift distributions used in the model are 
co rrect. We h ave c ompared our models with measurements 
of lVieroetaTI (1201 2h between SPIRE bands (see Fig.|7]i. All 
versions of our model agree with the data, except a systematic 
2cr tension for the point at { - 700 in the three bands which 
could be due to cosmic variance (note that one 2 o- outlier for 
18 poi nts is statistically expected). The trends of iHaiian et al.l 
(1201 lb (see Fig.[8l) are also well reproduced, but our model is 
systematically lower than the data at ^ = 1500. However, the 
stability of BLAST at large scales is not as good as for SPIRE. 
A future analysis of the cross-correlation between SPIRE and 
ACT and/or between Planck bands will thus be useful to further 
investigate this discrepency. 



6.3. Cross correlation between CIB and CMB lensing 

In Fig.|9] we compare the predictions of the cross-correlation 
between CIB and CMB lensing with the measurements from 
iPlanck Collaboration et aP (l2013h . At £ >1000, all versions of 
the model agree well with the data. At { <1000, the various 
versions of the model tends to be be systematically higher 
than the data. Models B and C (with modified evolution at 
high redshift) are closer from the data, but over-predict the 
cross-correlation by Icr at 857 GHz and 545 GHz and 2cr at 
353 GHz. Future works will investigate this tension. 



at smaller scales at 250jum and 350/^m, where systematics could 
be due to background subtracti on. Ou r model disagrees with the 
measurements of Cooray et alj ( |2010|) . which are systematically 
higher More reliable measurements of the correlation function 
of SPIRE sources, controlling the systematic effects, are thus 
necessary to check accurately the validity of our model. 



7. Successes and limitations of our model 

In the previous sections, we presented a natural way of ex- 
tending the B12 model of infrared galaxies by linking their 
properties to their host halo. We tested the validity of this 
approach comparing the prediction of this extended model to 
the measured spatial distribution of both individually-detected 
infrared galaxies (where the measurements are reliable) and 
the unresolved background. This comparison shows the good 
predictive power of our approach, suggesting that our assump- 
tions are fair Using this model to interpret both the origin of 
CIB anisotropics (Sect.[8]l and the link between star-forming 
galaxies and dark matter halos (Sect. 19. 3] ) is thus legitimate. For 
simplicity, we will perform this analyses only with model C, 
which provides the best agreement with the data. 

However, some limitations of this model have to be kept 
in mind. Firstly, the SED of infrared galaxies at high redshift 
(z > 2) was not measured and it is not clear if the dust temper- 
ature will increase or decrease at high redshift, because of the 
uncertainties on the mass-metallic i ty and mass-attenuation rela- 
tion as discussed in iMagdis et al.l (l2012h . The good agreement 
between the model and the CIB anisotropics at long wavelengths 
suggests that our hypothesis of non-evolution is reasonable. The 
SED of low mass galaxies (M* < 10'" Mq) is not well con- 
strained, but the contribution of these low mass galaxies to the 
CIB is small and this scenario is thus hard to test. In addition, the 
specific star formation rate and the SMF of star-forming galaxies 
at very high redshift [z > 3) are also uncertain, but the effects of 
z > 3 galaxies on the CIB are pretty small and these data thus 
cannot accurately constrain the evolution of infrared galaxies at 
very high redshift. Finally, we assumed an universal relation be- 
tween stellar mass and UV attenuation, but this relation could 
also break at very high redshift, where the amount of metals is 
smaller and galaxies could thus contain less dust. For all these 
reasons, the predictions of the models at z > 2 must be inter- 
preted with caution. 



6.4. Clustering of resolved galaxies 

The clustering of bright resolved infrared sources is also an 
important test of our model. Fig.[TO] shows the comparison 
between the clustering measurements of various flux-selected 
samples and our model predictions. At 100 and 160;um, our 
model agrees well with the HerschelfPACS measurements of 
iMa^Uocchetti et al.l (l201lh . In i/ersc/je//SPIRE bands, the com- 
parison is more tricky. There are str o ng dis agree ments between 
the m easurements of ICoorav et al.l (1201 Ol) and iMaddox et al.l 
(I2OIOI) . In fact, measureming the correlation function of SPIRE 
sources (250, 350, and 500 //m) is v ery difficult because these 
data are strongly confusion-limited dNguven et al.ll2010h . The 
background is thus hard to estimate and the completeness of the 
catal og can vary depending on the local source density. We agree 
with lMaddox et alT(l2010l) at scales larger than 0.1 deg, but not 



8. Where do the CIB and its anisotropies come 
from? 

In this section, we discuss the predictions of the model concern- 
ing the redshift distribution of the various signals we studied 
(CIB auto-spectra, CIB cross-spectra, CIBxCMB lensing). We 
base our analysis on the model C as motivated in Sect.|2] 

8.1. Redshift distribution of CIB mean level and anisotropies 

The CIB SED can be predicted directly from the lBethermin et al] 
(l2012ah counts model. Fig.[TT] (left upper panel) shows the 
spectral energy distribution of the CIB and its decom position 
per redshift slice. Compared to iBethermin et al] (1201 ll) model, 
this new model has a slightly higher contribution of 1 < z < 2 
sources. This new model agrees well with both absolute 
measurements and the total contribution of infrared galaxies 
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Fig. 11. Predicted redshift distribution of the CIB and its anisotropics. Upper le ft panel: CIB SE P and contribution per redshift shce 

(colo r ed sohd lines), black squares: total extrapo lated CIB from deep number counts jTeplitz et alj|201 1 : Bethermin et al. 2010a; Bertaet alj 
1201 iLlBethermin et al.ll2012d : lzemcov et al.ll20Ta) . Cyan solid line: Absolute CIB spectrum measured by COfiE/FIRAS (La gache et al...200a') . 
Green triangles: absolute CIB meas urements performed by COSii/DIRBE at 100 //m, 140 yum, and 240 yum (updated in Dole e t aLll2006h . Yellow 
diamond: absolute measu rements of [Penin et al.] l l2012ah at 160 yum with 5p//zer/MIPS. Orange arrows: upper limits derived from opacity of the 
Universe to TeV photons iMazin &R aue 2007). Lower left panel: Normalized redshift distribution of the mean level (up), Poisson anisotropics 
(middle), large scale anisotropics at ^=1000 (bottom). Various colors correspond to various bands. Right panel: Contributio n of v arious redshifts 
to the CIB (solid line) an d comparison with lower limits derived by stacking fromBerta et al. (2011) at lOOyum, Jauzac et al.: d20Ill) at 160 ^m and 
iBethermin et alj ( 1201 2ch at 250 yum, 350 yum, and 500 /jm (arrows). The dashed lines are the model predictions taking into account the selection 
used to derive the lower limits. 



extrapolated from the number counts. 

The three lower left panels of Fig.[TT] show the redshift 
distribution of the intensity of the CIB (top), but also of its 
anisotropics at small (middle) and large scales (bottom). In order 



to allow an easier comparison between bands, we normalized 
the redshift distributions in order to have J d(vly)/dzdz - 1 
and JdC(/dzdz = 1. Between 100 jum (3000GHz) and 850//m 
(353 GHz), the redshift distribution evolves strongly toward 
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Fig. 12. Predicted redshift distribution of auto and cross power spectra at large ({ = 1000, left) and small scale (Poisson level). We focus only 
on three frequencies: 857 GHz (350//m), 545 GHz (550/jm), and 217 GHz (1382 /jm) to illustrate the difference between cross-spectrum between 
two distant wavelengths and the auto-spectrum at an intermediate wavelength. 

Table 1. Predicted correlation of GIB anisotropies between bands at small scale (Poisson level, upper part) and large scale (1=1000, lower part). 



Correlation between bands at small scales (Poisson) 
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1.000 
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0.333 
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0.940 
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1.000 


Correlation between bands at large scale (1=1000) 
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217 GHz 
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0.997 


1.000 






143 GHz 


0.538 


0.904 


0.975 


0.996 


0.999 


1.000 




100 GHz 


0.575 


0.919 


0.981 


0.996 


0.997 


0.999 


1.000 



higher redshift. At longer wavelengths, there is almost no 
evolution, because all the sources below z ~ 1 are seen in the 
Rayleigh- Jeans regime and the color is roughly the same at all 
redshifts. This trend is seen for both intensity and anisotropies, 
regardless of scale. There are, however, small differences 
between these three quantities. For instance, the redshift 
distribution of the Poisson leveO at lOO/^m is dominated by 
Z < 0.5 sources, when a significant fraction of f = 1000 
anisotropies is caused by galaxies at z~l. This is due to the fact 
that the Poisson anisotropies are dominated by a small number 
of low-z bright galaxies just below the flux cut. Anisotropies 
at large scale are dominated by normal star-forming galaxies 
at z~l which dominate the background and the 2-halo term of 
anisotropies. 

Finally, we checked that the redshift distributions predicted 
by our model agree with the lower limits derived by stacking 
of 24 ^m sources . We co mpare with the low er limits from 
lBertaetalJ ( l2011h at 100/um. |jauzac et all ( 1201 Ih at 160 //m and 



We used the same flux cuts as in Sect.|6] 



iBethermin et al.1 ( 1201 2d ) at 250 //m, 350 /^m, and 500 //m, which 
well agree with our model. We also simulate with our model 
the selection corresponding to these various works and found a 
good overall agreement below z~2.5. At larger redshift and for 
a 24yum selection, these predictions are underestimated. This 
is expected because our SED templates take only into account 
the dust emission, when /} < 8 fim rest-frame emission are sig- 
nificantly affected by stellar emissions. We thus underestimate 
the number of 24yum detections and consequently also their 
contribution to the CIB. 



8.2. Redshift distribution of CIB cross-correlation between 
bands 

Cross power spectra between bands provide complementary 
information to the auto-spectra. Fig. [12] shows the redshift 
distribution of a selection of cross and auto-spectra. At both 
large {£ = 1000, left) and small scale (Poisson), the redshift 
distribution of the cross-spectra between two distant bands (e.g. 
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In contrast, the correlation (CJ"'/ JCJ" x Q'"') between 

bands provides a useful test of the validity of the models 
(especially their SEDs and redshift distributions). This corre- 
lation varies with scale. We thus focused on Poisson level and 
E - 1000. The results are summarized in Table[T] The level of 
correlation is also important for cosmic microwave background 
(CMB) experiments to know wether or not the CIB at high 
frequency is a good proxy of the CIB emission at low frequency. 
For instance, the correlation between 857 GHz and 143 GHz is 
0.904 at f = 1000, but only 0.75 for the Poisson level. This 
correlation is predicted to rise to 0.996 between 353 GHz and 
143 GHz (0.963 for the Poisson level). According to our model, 
the CIB at 353 GHz thus provides a very good proxy for CIB in 
the CMB channels. 



8.3. Redshift distribution of CIBxCMB tensing signal 

The cross-correlation between CIB and lensing provides supple- 
mentary information probing a slightly different redshift range. 
Fig. [13] shows a comparison between the predicted redshift 
distributions of the CIB auto-spectrum and the cross-spectrum 
between CIB and CMB lensing for various bands. While the 
redshift distribution of the auto-spectrum evolves strongly 
between 100 (3000GHz) and 850yum (353 GHz), this evolution 
is weaker for the cross-correlation with CMB lensing which 
mainly probes the z - I - 2 redshift range. CIBxCMB lensing 
is thus a good probe of the mean SED of galaxies contributing 
to the CIB at this epoch. However, apart from lOOyum band, the 
redshift range probed by CIBxCMB lensing is relatively similar 
to that of the CIB power spectrum, but it is an independent 
and more direct probe of the link between dark matter and 
star-forming objects, because the lensing of the CMB by the 
large scale structures is a very well known and modeled physical 
process. 



9. Star formation and dark matter halos 

In this section, we will discuss the prediction of our model in 
the context of our current understanding of galaxy formation. 
We use model C, as discussed in Sect.|2] 



Fig. 13. Comparison between the predicted redshift distributions of 
CIBxCIB and CIBxCMB lensing signal in various IRAS/Planck bands. 
Red solid line: Redshift distribution of the cross-correlation between 
CIB and CMB lensing at 1=1000 in various bands. Blue dashed line: 
Redshift distribution of CIB auto-power-spectrum. 



857 and 217 GHz here) is pretty close to the redshift distribution 
for an intermediate band (e.g. 545 GHz here), except for a 
difference of Poisson level at very low z caused by different flux 
cuts. Considering the quite fine spectral coverage of our data, 
we thus cannot claim that cross-spectra probe different redshift 
ranges. 



9.1. Evolution of infrared-light-to-mass ratio with redshift 

The infrared-light-to-halo-mass ratio is a key ingredient of our 
CIB model. Fig. fT4l shows this relation for various hypotheses. 
We consider two different definitions of the infrared luminosity 
inside the halos: the infrared luminosity of a central galaxy ly- 
ing exactly at the core of the main-sequence (top panels) and 
the mean infrared luminosity of central galaxies (lower panels) 
hosted by halos of a given mass. This second quantity takes into 
account the fact that a fraction of galaxies are passive, while the 
first quantity takes into account only the star-forming galaxies. 
We also use both the instantaneous halo mass (left panels) and 
the halo mass at z=0. The conversion between instantaneous 
halo mass and halo mass a t z=0 is performed ass uming an ac- 
cretion following the fits of iFakhouri & Mai (1201 Ol) to their own 
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Fig. 14. Relation between infrared luminosity of galaxies and their halo mass. The top panels correspond to the infrared luminosity for a central 
galaxy exactly on the main-sequence. The bottom panels shows the mean infrared luminosity of the central galaxies. The left panels shows results 
as the instantaneous halo mass and right panel as a function of halo mass at z=0. The dashed lines correspond to the limit between normal galaxies 
(i.e. < 10"Lo), LIRGs(10" < Ljr < IO'^Lq), ULIRGs (10'' < Ljr < 10'^^ Lq) andHyLIRGs (> 10'^ Lq). 



numerical simulations of dark matter; 

,1.1 



<M,,)=46.1xMo/yr 



Ml, 



x(l + l.lk) 






(34) 



The mass at z=0 (Mi,{z) + [,,_, {Mi,)dt ) is a convenient 
quantity to use when discussing the evolution of structures with 
time. 



The shape of the infrared luminosity as a function of 
instantaneous halo mass evolves little with redshift. However, 
its normalization varies strongly. For the same halo mass, 
infrared luminosity is much more intense at high redshift. 
This results will be interpreted further in Sect. 19. 31 At low 
mass, the infrared luminosity decreases quickly for several 
reasons. The main cau se is the low star formation efficiency 
in low -mass halos (e.g. iBehroozi et al]l2012al : iLeauthaud et al.l 
I2OI2I) . In addition, low-mass galaxies are less attenuated than 
massive ones, such that only a small fraction of the UV photons 
emitted by young s t ar is reprocessed in to infrared light (e.g. 
iPannella et all 120091; ISobral et all l2013h . At high mass, we 
found a sub-linear relation, even if we consider only objects on 



the main-sequence. This implies that mechanisms of feedback 
are pretty efficient in massive halos. This trend is stronger if we 
consider the mean infrared luminosity taking into account the 
quenched galaxies. At z<l, the infrared luminosity decreases 
when halo mass becomes larger than IO'^'Mq, and our model 
thus predicts a very weak star-formation in dense environmen t 
at low redshift in agreement with e.g. iFeruglio et alj ( l2010l) . 
The contribution of satellites is small except in massive halos 
(> 10'^ Mq), where the central is inefficient in forming stars. 

The infrared luminosity as a function of halo mass at z=0 
exhibits a strong downsizing. The progenitors of 10'"* M© 
initially host a very intense star formation rate (ULIRG regime, 
10'~ < LiR < lO'"* L0) at high redshift (z=5) which continues 
to grow up to z~2. Then, the infrared luminosity decreases 
because star formation is less efficient in massive halos (see 
Sect. |9.3t . The less massive halos need more time to ignite 
intense star formation. Dark matter halos similar to that of 
MilkyWay (~10'~ M© at z=0) host significant infrared emission 
only below z= 1 . 

We can compare the predictions of our model with es- 
timates based on clustering of various samples of galaxies. 
Sub-millimeter galaxies (defined here to have 5 850 > 3mJy) 
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Fig. 15. Contribution of various halo mass to the star formation history. Upper panels shows the contribution of each shce to star formation 
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right panels. 



are essentially ULIRGs at z~2.4 (e.g. IChapman eFaH l2003h 
and the B12 model pre dicts that th e se obj ects generally lie 
on the main-sequence. iBlain et aP (l2004l) showed that the 
correlation length of SMGs is compatible with typical host halos 
of 10''' Mq, which is consistent with our model prediction of an 
infrared luminosity of 4 x lO'^L© for a main-sequence galaxy 
in a lO'^ Mq halo at z=2.4 (see Fig.fl4li. We can also check the 
consistency of our resul ts with measuremen ts based on Spitzer 
24-jjm observations. Magliocc hetti et all (l2008h measured 
the clustering of two samples of sources with ^24 > 400;uJy 
with a mean redshift of 0.79 and 2.02, respectively. This flux 
selection corresponds, at the mean redshift of each sample, 
to an infrared luminosity of 2.6 x lO" Lq and 2.6 x lO'^L©, 
respectively for a main-sequence SED. For an object exactly 
on the main sequence, this corresponds to an instantaneous 
halo mass of 2.7 x IO'^Mq and 1. 2 x lO'^Mp, in agr e emen t 
with the minimal mass found by iMagliocchetti et all (l2008h 
of 0.8+g;^ X IO'^ Mq and 0.6+°-^ x IO'^Mq, respectively. 
iFarrah et al.l (l2006l) measured the halo mass of 5.8 yum IRAC 
peakers with 24 //m detections and found a Icr range for log(M/,) 
of 13.7-14.1. This population of massive star-forming galaxies 
has a mean redshi ft of 2.017 and a mean infrared luminosity of 
8.9±O.6L0 (Fiole t et alj|2010t) . which is associated in our model 
tolog(M/,)=14.1. 



9.2. Contribution of various iialo masses to star formation 
liistory 

The total star-formation rate of both central and satellite 
galaxies hosted by a given halo can be derived from our model. 
In combination with the halo mass function, we can then 
compute the contribution of each halo mass to the star formation 
rate. The results of this analysis are presented Fig. [15] If we 
consider instantaneous halo mass (left panels), the bulk of 
the star formation is hosted by halos between 10"^ Mq and 
IO'-^'^Mq regardless the redshift, suggesting the existence of a 
characteristic mass which favors star-formation. The fractional 
contribution of each mass sUce is inferred to have evolved 
slowly with redshift, except at z > 3. However, there are 
significant uncertainties in the behavior of the galaxies at high 
redshift and these results should be taken with caution. 

Because of the growth of dark matter halos, the same 
halo mass at various redshifts corresponds to progenitors of 
different type of halos. In order to take this into account, 
we also considered slices of z=0 halo mass (right panels). 
The star formation rate density is predicted to have been 
successively dominated by progenitors of massive clusters 
{Ml, > lO'"*'^ Mq) at z > 4.7, small clusters and large group 
(10'^-^ < M,, < lO'^^-^Mo) at 2.8 < z < 4.7, small groups 
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Fig. 16. Instantaneous star formation efficiency (defined here as SFR/BAR, where BAR is the baryonic accretion rate) as a function of halo 
mass at various redshift predicted by our model. The predictions from model C are plotted with a solid line. At z = 5 model A predicts a different 
efficiency distribution and is represented by a dashed line. Upper panels shows this efficiency only for main-sequence galaxies and lower panels 
for mean efficiency. We use slice of instantaneous halo mass in the left panels and mass at z=0 in the right panels. 



(10'--^ < M/, < W^-^ Mo) at 0.5 < z < 2.8, and Milky Way-like 
halos (10"-5 < Mh < 10'" Mo) at z < 0.5. Star-formation 
thus initially occurs in the progenitor of the most massive 
halos before becoming less efficient, and then propagating to 
less massive halos. This strong downsizing explains why the 
star formation rate densities at redshift 1 and 2 are similar. 
z > 2 infrared luminosity density is dominated by a small 
density (~ 3 x lO^^Mpc '*) of ULIRGs and z ~ 1 is domi- 
nated by a larger den sity (~ 3 x lO^^Mpc"^) of LIRGs (e.g. 
iBethermin et al.ll201 lb . 



9.3. Efficiency of star formation as a function of inaio mass 
and redsliift 

To conclude, we derive an estimate of the instantaneous star 
formation efficiency (ISFE) in different halos by computing the 
ratio between SFR and baryonic accretion rate (BAR). Here, 
we assume the BAR to be the total matter accretion as given 
bylFakhouri & Ma (2010) multiplied by the universal baryonic 
faction (BAR = {M,,} x 0/,/Q,„ with <M/,) defined in Eq.l34ll. 
This simplified definition neglects the fact that the gas is not 



transformed instant aneously into stars and large gas reser voirs 
are pres ent in (e.g. Daddi etal] 120101: iTacconi etal] 1201 2h and 
around (ICantalupo et al.ll2012h high-z galaxies. Fig. [161 shows 
predictions for the variation of this efficiency as a function 
of halo mass at various redshift. Star-formation efficiency in 
main-sequence galaxies as a function of instantaneous halo mass 
(upper left panel) slowly evolves between z=0 and z=2 with a 
slight increase of the mass of maximum efficiency with redshift 
(from 8 X 10" Mo at z=0 to 3 x 10'^ Mo at z=2). At z=5, this 
maximal efficiency is smaller and the halo mass where it occurs 
is higher. The increase of the mass of maximum efficiency 
could be due to a delay in the ignition of star-formation activity, 
because the BAR in ~10'~Mo halos is higher at high redshift. 
However, large uncertainties exist at high redshift where we 
have few constraints on star-forming galaxies. For instance, if 
we use model A (flat sSFR at z > 2.5) instead of model C (rising 
sSFR) to perform our calculations this efl'ect is much smaller 
(dashed line in Fig.fTSIl. Below the mass-scale at which peak 
efficiency is reached, the gravitational potential of the halo is 
lower and supernova feedback is pr obably sufficiently strong 
to rem ove gas from the galaxy (e.g. ISilkl 120031 : iBertone et al.l 
l2005h . At higher masses, the slow decrease of the star formation 
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efficiency could be caused by the transition from cold streams 
below M/, ~ IO'^M q to isotropic cooling above this mass (e.g. 
Ibekel & Birnboiml F2006; Faucher-Giguere et al. 2011). At high 
mass, the cooling time of gas becomes much longer than the 
free fall t ime because of the hot atmosphere in the massive 
halos (e.g. lKeres et alj|2005l: iBirnboim et alj|2007l) . Contrary to 
the ISFE in main sequence galaxies, the mean ISFE (lower left 
panel), that takes into account the fact that a fraction of galaxies 
are quenched, exhibits a strong break at M/, > 10'^ Mq at z = 
and a moderate break at the same mass at z = 1. This break 
could be caused by a suppression of isotropic gas cooling by 
energy injection in the halo atmosphere by ac ti ve galactic nuclei 
(AGN) activity (eg. iCattaneo et al.i 120061: ISomerville et alJ 
l2008t IOstrikeretani20M 

Regardless of the redshift and halo mass, the SFR in main- 
sequence galaxies is always lower t han the BAR ont o the dark 
matter halos (see Fig.O. For a IChabrieil ( l2003h IMF, the 
maxi mal effic i ency i s ~0.7. AppendixICl discusses the case 
of a ISalpeteiJ (Il955h IMF. This is consistent with the stan- 
dard assumption that main-sequence galaxies host secular star- 
formation. This is not the case for episodic starbursts, which on 
average forms 4 times more stars for the same halo mass than 
main-sequence galaxies and thus transform much more gas into 
stars than they receive from cosmic accretion. They thus tend to 
exhaust their gas reservoir rapidly. The fact that main-sequence 
galaxies do not need to store gas to fuel their star-formation 
could seem in contradiction with the fact that the specific halo 
growth rate (sHGR) is larger b y a factor of two t han th e sSFR 
around z = 2, as pointed out bv 'Wein mann et all ( |20I1|) . How- 
ever, the ratio between sSFR and sHGR is: 



(35) 



SFR/M^ ^ r^MblM,, ^ Ah 
Ml, /Ml, M^/Mi, ^M^' 

where rj is the instantaneous star formation efficiency (SFR/M/,), 
Ml, the baryonic accretion into the halo. The ratio between 
sSFR and sHGR can thus be much higher than unity in halos of 
~ IO'^Mq, because, even if rj is lower than 1, the ratio between 
baryo nic mass and stellar mass is ~4 (e.g. Leauthaud et al. 
l20I2h . because of the low efficiency of conversion of baryons 
into stars in the past (when the halos had a lower mass and thus 
a lower ISFE). 

We can also discuss how the efficiency varies as a function 
of z=0 halo mass (right panels). We predict a strong evolution 
of the typical halos where the ISFE is maximal. The progenitors 
of massive halos form stars very efficiently at high redshift, 
while MiUcy-Way-like halos are very inefficient. The opposite 
trend is expected at low redshift. This picture is consistent with 
the strong downsizing of the star-forming galaxy population 
discussed in t h e previ ous sections, but also with the work 
iBehroozi et al.l (1201 2bh . Their analysis was based only on 
the evolution of the SMF and SFR were derived assuming a 
single galaxy population. Our approach is based on the infrared 
observations which directly probe the SFR in galaxies and 
which take into account the diversity of galaxies using three 
distinct populations: secularly star-forming galaxies on the 
main-sequence, episodic, merger- driven starburst s and pa ssive 
elliptical galaxies. Compared to IBehroozi et al.l (l20I2ah . our 
mean ISFE is much lower in local massive halos, which almost 
exclusively host quenched galaxies, but we agree with their 
estimate if we take only a main-sequence population (more 
consistent with their single galaxy population). 



Normal spiral galaxies at z = 0, LIRGs at z = 1 , and ULIRGs 
at z = 2 dominate the star-formation density at these redshift and 
are essentially main sequence galaxies (Sargent et al. 2012). All 
of these are hosted by halos of similar mass associated, which are 
characterized by a very efficient conversion of accreted baryons 
into stars. The huge difference between their star formation rate 
can thus be explained by the accretion, which is stronger at high 
redshift. These objects can thus be viewed as different facets of 
the same universal process of secular star formation. 

10. Conclusion 

We have studied the connection between star formation and dark 
matter halos focusing on infrared observations. We developed a 
new modeling approach based on the 2SFM framework, which 
was already able to reproduce su ccessfully the infrar ed luminos- 
ity function and number counts dSargent et al.ll20I2l B12). This 
framework links stellar mass with star formation and infrared 
properties, assuming two different modes of star formation in 
secularly star-forming galaxies and episodic starbursts. We ex- 
tended this formalism to the connection between stellar mass and 
halo mass using the technique of abundance matching. Our for- 
malism account for the fact that a significant fraction of massive 
galaxies are passive and do not form stars (especially at low red- 
shift). 

- We developed a method to compute the CIB anisotropics 
(including power-spectra between different frequencies), the 
cross-correlation between CIB and CMB lensing, and auto- 
correlation functions of bright resolved galaxies using the 
prescription of the 2SFM formalism. To perform this compu- 
tation, we produced effective SEDs of all galaxies at a given 
redshift which are available onlinq^- 

- We find that a slowly rising sSFR at z > 2.5 and a quenching 
of satellite of massive passive galaxies at low redshift (our 
model C), matches reasonably well the infrared data suggest- 
ing that this model is a valid description of the link between 
infrared galaxies and dark matter halos. The other versions 
of the model with flat sSFR and no quenching of satellites 
show some little discrepancies with the power-spectra, but 
cannot be conclusively ruled out. 

- Our model is able to predict the redshift distribution of CIB 
anisotropics. We found that the mean redshift where the CIB 
is emitted varies strongly between lOOyum and 850 yum but 
only little at longer wavelengths. Consequently, the CIB 
anisotropics in the various bands above 850 yum are strongly 
correlated (> 0.9). 

- We found a quick rise of the far-infrared-light-to-mass ratio 
with redshift in M/, > 10'^ Mq halos and a strong break at 
lower mass at all redshifts. We found that more than 90% of 
the star formation is hosted by halos with masses between 
lO"'^ and 10'^'^ Mq at all redshifts. The progenitors of 
clusters {Mi,{z = 0) > 10'^-^ Mq) host the bulk of the star- 
formation at z >3. Star formation activity then propagates 
to groups (10'--^ < Mi,{z = 0) < 10'^-^ Mg) at 0.5 < z < 3 
and Milky Way-like halos (10"-^ < Mi,(z = 0) < IO'^-^Mq) 
at z < 0.5. We also found that there is a characteristic halo 
mass (~ 10'^ M0) where the star formation efficiency is 
maximal (~70% at all redshift). The large difference of 
SFR in galaxies dominating the background at low and high 
redshift would thus be driven by a difference of accretion 
rate in the halos close to this mass. 
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Our simple modeling framework is very efficient in explain- 
ing the current observations of the infrared Universe. However, 
future large sub-millimeter surveys (e.g. NIKA, CCAT) will 
resolve the bulk of the CIB into individual sources and will 
probably improve measurements of the clustering properties of 
infrared galaxies. We expect that deviation from our model will 
appear at small scales, where environmental effects could have 
a strong impact, thereby revealing a more complex and varied 
infrared Universe. 
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Appendix A: Conversions tables 

Appendix A. 1 : Wavenumber, multipole, and angle 

On small scales {Spitzer, Herschel), where curvature of the sky is 
negligible, people uses in general wavenumber (k) in their mea- 
surements of CIB power-spectrum. This is not the case for large- 
scale measurements (SPT, ACT, Planck) for which people used 
mutlitpole {(). The conversion between the two is just C - Ink. 
This also corresponds to a characteristic angular scale - n/t. 
Table lA. ll provides conversion for the range of values used in this 
paper 
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Table A.l. Conversion between wavenumber, multipole, and angle 



Multipole 


Wavenumber 


Angle 


{ 


k 


e 




arcmin"' 


degree 


1 


0.00005 


180.00000 


2 


0.00009 


90.00000 


5 


0.00023 


36.00000 


10 


0.00046 


18.00000 


20 


0.00093 


9.00000 


50 


0.00231 


3.60000 


100 


0.00463 


1.80000 


200 


0.00926 


0.90000 


500 


0.02315 


0.36000 


1000 


0.04630 


0.18000 


2000 


0.09259 


0.09000 


5000 


0.23148 


0.03600 


10000 


0.46296 


0.01800 


20000 


0.92593 


0.00900 


50000 


2.31481 


0.00360 


100000 


4.62963 


0.00180 



Table A.2. Conversion between wavelength and frequency for various 
passbands used in this paper 



Wavelength 


Frequency 


/im 


GHz 


24 


12500 


100 


3000 


160 


1875 


250 


1200 


350 


857 


500 


600 


550 


545 


850 


353 


1100 


272 


1382 


217 


2097 


143 


3000 


100 



(B.2) 

Since we assume Poisson statistics, the variance on the num- 
ber of galaxies in this bin equals n^^^"^®. For the co-variance 
between the fluxes at the two frequencies caused by this sub- 
population we have: 



of 



SorSB 



— L> v^^ v' ^ f^ijk — ^IR i'^ 



MS or SB ^ MS 01- SB 



X'Vjk, (B.3) 



because the fluxes in the two bands are perfectly correlated for 
sources with the same SED. Finally, we sum over the entire pop- 
ulation (all LiR, z, and (L'^)) to compute the level of the Poisson 
term: 



cz.- Z ZZZ< 

(MS,SB|LiR,i Zj {U)t 



SorSB 



which in integral limit becomes: 



<'= ff Z r /'MSo.SB«t/)k) P" 



(B.4) 



"«£/>,.') 



(MS.SB) 

d A^MSoi-SB ,2 MSorSB//Tr\ ,^ ^ MSotSBz/ttx n .r J/Tr\^^ 
clLijiu V 

(B.5) 



Appendix C: Star formation efficiency in the case of 
a Salpeter IIVIF 



In Sect. 19. 31 we computed the SFE assuming a IChabri er ('2003J) 
IMF. T he results are slightly different if we assume a .Salpeter! 
(1 19551) IMF (see Fig. lB.lb . In this case, the galaxies close to 
the mass of maximum ISFE form more stars than they accrete 
baryons and large gas reservoirs are required to allow the 
secular star-formation in these objects. These reservoirs could 
have been replenished during the phase of low star-formation 
efficiency, when the halo was less massive. 



Appendix A. 2: From wavelengths to frequencies 

Infrared astronomers use wavelengths in yum, the CMB commu- 
nity frequencies in GHz. For quick reference, we here provide in 
Table lA.2l the conversion between these two conventions for the 
passbands discussed in this paper. 



Appendix B: Computation of the cross-Poisson 
term 

The number of sources per steradian ritj with an infrared lu- 
minosity in the interval [LiR,i,LiR,i + ALiri] and a redshift in 

[Zj,Zj + Azj]is: 



dV 



d^N 



dz dLiRdV 



X 



X ALiR.iAzj. 



(B.l) 



We also need to separate the galaxies by mode of star formation 
(main-sequence or starburst) and (U) parameter. The number of 
main-sequence or starburst galaxies n^^°'^^ in the same Lir and 
z bins and with (U) in [([/)*, {U}k + A{U}k] is: 



^MSorSB 
^ijk 



dV 



X 



d^N 



dz dL\^dV 



X ALiR,iAz; X PMSorSB«t/» X A{U)k. 
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Fig. B.l. Same figure as[T6]but assuming a lSalpeten il955f) IMF. Tlie black daslied line corresponds to a SFE of 1. 
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